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We develop a method to compute the fermion contribution to the vacuum polarization 
energy of string-like configurations in a non-abelian gauge theory. This calculation has 
been hampered previously by a number of technical obstacles. We use gauge invariance of 
the energy and separation of length scales in the energy density to overcome these obstacles, 
^vq , We present a proof-of-principle investigation that shows that this energy is small in the MS 

renormalization scheme. The generalization to other schemes is straightforward. 
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INTRODUCTION 



Various field theories suggest the existence of string-like configurations, which are the particle 
physics analogues of vortices or magnetic flux tubes in condensed matter physics. Often they 
are called cosmic strings to distinguish them from the fundamental variables in string theory 
and to indicate that they stretch over cosmic length scales. If they exist, cosmic strings, can 
4h ■ potentially have significant cosmological effects. We refer to ref. [l| for a recent review on the 

physical implications of strings in the standard model and beyond. 
CN ■ In the standard model of particle physics, string solutions are not topologically stable 

and thus can only be stabilized dynamically. In exploring the existence of cosmic strings, it is 
therefore important to be able to accurately calculate their energies. Here we will develop and 
apply a method to compute the fermion contribution to the leading order quantum correction to 
the energy, the so-called vacuum polarization or Casimir energy. In a large Nq scenario with many 
internal degrees of freedom, the fermion contribution dominates that of the bosons. 
' A number of previous studies have investigated quantum properties of string and vortex con- 

figurations, invoking either approximations, simplified configurations, or lower dimensions to cope 
with technical difficulties. Naculich jsj has shown that in the limit of weak coupling, fermion 
fluctuations destabilize the non-abelian Z-string. The quantum properties of Z-strings have also 
been connected to non-perturbative anomalies 0]. The fermionic vacuum polarization energies 
of QED flux-tubes, and abelian flux-tubes more generally, were investigated using heat-kernel 
methods 0, 0], world line numerics as well as the phase shift method 10]. The heat-kernel 



method was also used to study self-dual vortices [HI]. This method and the world line approach 
limit renormalization to the subtraction of the divergences in the heat -kernel expansion, thereby 
obscuring the connection to perturbative renormalization. As we will explain later, the phase shift 
approach is capable of making straightforward contact with any renormalization condition that is 
formulated in terms of (momentum space) Green's functions. In lower dimensions the ultra-violet 
divergences are less severe, which made the computation feasible for the case of two spatial dimen- 
sions [13], while the obstacles that arise for the physical case became soon obvious [13!]. A first 
attempt at a full calculation of the quantum corrections to the Z-string energy was carried out in 



ref. [TJ]. Those authors were only able to compare the energies of two string configurations, rather 
than comparing a single string configuration to the vacuum; these limitations arise from subtleties 
of the renormalization process that we address in this paper. Also, the contribution of bosonic fluc- 
tuations to the vacuum polarization energies have been estimated for vortex configurations using 



the heat-kernel method [15(] and string backgrounds within the phase shift approach [la ], which 
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we will use here for the fermionic fluctuations. Stability of cosmic string currents was considered 
in 



17] 



We begin by expressing the vacuum polarization energy as the renormalized sum over (half) the 
change of the single particle energies caused by the localized background. Previously, the fermionic 
vacuum polarization energy of strings has been computed for the case of two spatial dimensions [l~2| . 
The main purpose of the present paper is to extend this computation to the physical case of three 
spatial dimensions, yielding the vacuum polarization energy per unit length of the string. Though 
that extension seems straightforward since the string is translationally invariant with respect to this 
additional coordinate, a number of obstacles arise. They are mainly related to the more complex 
structure of ultra-violet divergences. 

It is well established that the vacuum polarization energies of extended background field config- 
urations, such as solitons or vortices, are unambiguously obtained from a momentum integral that 



involves the derivative of the phase shifts in the potential generated by the background field 18]. 
These phase shifts measure the distortion of the spectrum of quantum fluctuations caused by the 
background. If an object is translationally invariant with respect to a subset of the coordinates, we 
can use the phase shifts calculated in the nontrivial dimensions combined with appropriate kine- 



matic coefficients to describe the full spectrum of fluctuations [191 ] . These coefficients vary only 
with the number of trivial dimensions, but not with the background field. In addition to integrat- 
ing the result over the magnitude of momentum k, we also need to sum over angular momentum 
channels. It has previously been shown that for string-type configurations, these two operations 
are not absolutely convergent, and inappropriately exchanging them may yield an unphysical con- 
vergence 131 ] . In D = 2 + 1 spacetime dimensions, the problem is mitigated because the relevant 
Feynman diagram is manifestly finite. In D = 3 + 1, the complicated structure of the divergent 
third and fourth order Feynman diagrams makes the exchange unavoidable. This problem is not 
specific to the string, and can only be avoided by analytic continuation to complex momenta. This 
procedure requires a careful construction of the Jost function for the Dirac scattering problem. On 
the other hand, the integration over imaginary momenta automatically includes the bound state 



contribution to the vacuum polarization energy [18|, 12 



A second, more severe obstacle is specific to the string configuration. The method for com- 
puting the vacuum polarization energy also requires the evaluation of Feynman diagrams whose 
external legs are given in terms of Fourier transforms of the background field. In this calculation, 
the leading terms of the Born series are subtracted and the corresponding contributions to the 
vacuum polarization energy are added back in as (renormalized) Feynman diagrams. We stress 
that the energy is not computed from this perturbation expansion; rather, the advantage of our 
approach is that this expansion is used to make contact with standard procedures of perturbative 
renormalization, which in particular allows the implementation of any renormalization scheme that 
is formulated for the momentum space Green's functions of the theory. The problem that arises for 
the string is that the configuration does not vanish at spatial infinity, and so its Fourier transform is 
ill-defined. While gauge-invariant combinations of the Higgs and vector fields do vanish at spatial 
infinity, the Feynman series is not gauge invariant term by term. To circumvent this problem, we 
need to introduce a return string, localized at a some distance po that is large compared to the 
typical extension w of the physical string. This return string unwinds the physical string so that 
the Fourier transforms and thus the individual Feynman diagrams give well-defined functionals. 
Even though the vacuum polarization energy E v3jC is a non-local functional, for pq ^> w there 
should be a separation of scales, 

-Evac = E ps (w) + E rs (p ) (1) 

for the energies of the physical string (E ps ) and that of the return string (E rs ), since the overlap 
between the associated densities vanishes. Then computing E vac and finding E TS (po) for a special 
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case, yields E ps (w), the quantity that we are seeking. A similar approach has been previously used 
to compute the vacuum polarization energy of electromagnetic flux tubes. In that case dimensional 
analysis indicated that linip _>. oo E TS (po) = 0, which was then confirmed by directly calculating the 
energy density [13]. Here the situation is more complicated, due to the presence of the Higgs 
field. A naive analysis indicates that in general E TS (po) actually grows quadratically with pq. If 
the Higgs field of the return string stays on the chiral circle (meaning that its magnitude is fixed 
at its vacuum expectation value, so that only its isospin orientation changes in space), E ts (pq) 
remains constant. However, there is no a priori reason for it to vanish. As an aside we note 
that, in contrast to the QED case, the classical energy of the return string does not vanish as 
Pq — > oo. So we need to find the energy of the return string, which then is subtracted from eq. ([I]) 
to separate the energy of the physical string. Unfortunately, the general return string background 
induces potentials in the scattering problem that behave like 1/p 2 at small distances. Though the 
corresponding scattering matrix can still be computed, the evaluation of the Jost function remains 
obscure because the bounds needed to prove the analytic structure from the iterative solution are 
no longer satisfied. One way out of the dilemma would be to compute the return string energy 



within a covariant expansion that is reliable for large po 2l|, but we have a better option. At least 
in a particular gauge (and we have to choose one anyhow to do the calculation) the small distance 
singularity disappears for return strings that have support only on the chiral circle. As a result 
we can calculate E rs in this case, and numerically we find that it vanishes as po — > co in the MS 
renormalization scheme. This is not in contradiction to the general analysis that yields a constant, 
because the non-zero term is omitted in MS. 

Another technical obstacle that did not emerge in the case of two space dimensions is the in- 
creased computational effort needed. In the case of three space dimensions, the third and fourth 
order Feynman diagrams induce logarithmic divergences. When folded with the Fourier trans- 
forms of the background potentials, these diagrams become higher-dimensional integrals that are 
cumbersome to treat numerically. We therefore introduce a fake boson field whose second order 



Feynman diagram possesses the identical logarithmic divergence [22], |23j . Since it is a solution to 
a well-defined scattering problem, we are ensured of the analytic properties of the resulting Jost 
function and we can deal with it on equal footing as the one from the Dirac problem. 

In this paper we have collected the ingredients for computing the vacuum polarization energy 
of the string. For simplicity we use the MS renormalization scheme and leave the use of physical 
(on-shell) conditions and dynamical questions such as the stabilization of the string by bound 
fermions to a forthcoming paper. 

To be specific, we consider a left-handed SU(2) gauge theory in which a fermion doublet \E' 
is coupled to the gauge field and the Higgs field cj). Since we are studying only fermionic 
fluctuations, we can ignore the self-interactions of the Higgs and gauge fields and simply assume 
the Higgs potential is minimized at a nonzero vacuum expectation value \(j)\ = v. We express the 

Higgs doublet 4> = ( ^ + ) as the matrix 



We then have the interaction Lagrangian 

c = m ltl D»p L y + m^d»p R y - / * Up r + &p l ) * , (3) 



where Pr,l = \ (1 ± 75) are projection operators on right- and left-handed components, respec- 
tively. The Higgs-fermion interaction is parameterized by the Yukawa coupling constant /, while 
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the gauge coupling constant is defined via the covariant derivative 



(4) 



Note that the vector fields are 2x2 matrices in (weak) isospace. To simplify the model, we have 
taken the Weinberg angle to be zero, so that there is no photon coupling, and assumed that the 
fermion doublet is degenerate in mass. 

Finally we list the counterterms that are required to compute the vacuum polarization energy, 



£ ct = Cl tr 7 [W^W^] + ^ tr, (D^ 



+ 2-trj 



$t$ 



+ 



tr/ 



(5) 



where the gauge field tensor is W, 



d^W u 



■dyW^ — ig [Wfj,, W u ] for the vector field. As indicated, 
the traces are with respect to isospin. Of course, the counterterms are just the typical gauge- 
invariant terms appearing in the classical Lagrangian for the Higgs and gauge fields. Note also 
that the C3 and C4 counterterms cancel quadratic ultra-violet divergences, while c\ and C2 cancel 
logarithmic divergences. In a model with only two spatial dimensions, c\ and C2 would be finite, 
though not necessarily zero 121 ] . 

The remainder of the paper is organized as follows. In the next section we introduce the string 
configuration and discuss the obstacles mentioned above in more detail. In section III we briefly ex- 
plain the phase shift method for computing the vacuum polarization energy. In particular we discuss 
the interface method for the case when the background configuration is translationally invariant 
in one coordinate. We also introduce the fake boson technique for simplifying the higher-order 
divergences. In section IV we explain the use of the return string on the chiral circle. Numerical 
results are presented in section V, and concluding remarks and an outlook for applications are 
given in section VI. 



II. THE Z STRING 

The string configuration only depends on the distance p from the symmetry axis (which we 
choose to be z) and the corresponding azimuthal angle (p. It is characterized by a non-vanishing 
angular dependence at spatial infinity. In matrix notation we write in temporal gauge (W° = 0) 

w ■ (t Jg{p) , ( sin(Ci) *cos(£i)e-^\ 

W = ^^—^{-icosi^ -sin(6) J 

JfliW \-"s(a) sin(6)e m V ' { ) 

This configuration is commonly called a Z-string, because the corresponding component Z oc 
tr i(Wt^) exhibits the spatial dependence of an abelian string. The radial functions fc(p) and 
fn(p) approach unity at spatial infinity while they vanish at the origin (p = 0). They are the 
typical profiles of the Nielson-Oleson string (24| . We parameterize these profile functions via 

—e- -(-e-Y 
f H (p) = 1 - e »fl and f G (p) = 1 - e \*>aJ . (7) 

Then the fields W and <3? are 0{p) as p — > and ambiguities resulting from an ill-defined azimuthal 
angle do not arise. We have also introduced a general winding number n for the string, though in 
the numerical calculation we will only consider n = 1. In ref. [121 ] we treated the angle < £1 < ^ as 
a variational parameter of the string configuration. We emphasize that £1 describes the orientation 
of the Higgs field on the chiral circle. It will prove to be very useful to introduce the return string 
by allowing it to be space dependent such that it vanishes at spatial infinity. 
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A. String in the Dirac Equation 



The Dirac equation that arises from the Lagrangian is dealt with in two steps. The dependence 
on the z-coordinate in which the string is invariant separates in a pure phase factor e~ %PzZ . In the 
second step we are thus left with a two-dimensional problem 1 



H 



+mf j 



H 





1 

-1 



-a • if a • ip 
a ■ (p —a ■ (p 



Ig 



+ s 



1 
-1 



Hp 



(8) 



where s = sin(£i) and c = cos(£i). The explicit matrices in the above equation refer to the spinor 
indices and <7j are the 2x2 Pauli matrices. Furthermore we have introduced the isospin matrices 



-ice 



riLip 



ice 



-imp 



and 



e m ^ 
e -irup o 



(9) 



If E is an eigenvalue of H, then the single particle energy is iy 7 E 2 + p^. Calculating the vacuum 
polarization energy requires us to integrate over p z in the framework of the interface formalism, 
which we review below. The sign degeneracy arises from the anti-commutator {if, 73} = 0. 
It is straightforward to extract the fermion-string interaction from eq. ([8]) 



Hmt — —fc 
2p 



-a ■ ip a ■ (p 
a ■ (p —a ■ (p 



Ig + m 



(cf H - 1) 



1 
-1 



+ sf H 



1 
-1 



Hp 



(10) 



In view of the asymptotic behavior of the radial functions /gh mentioned above, many of the 
problems with computing the vacuum polarization energy become immediately obvious. The in- 
teraction does not vanish at spatial infinity, but rather approaches a pure gauge configuration. 
This implies that there is no straightforward way to set up a Born series, whose necessity we will 
recognize in the following section. Furthermore we cannot compute any Fourier transformation 
of the interaction, which makes it impossible to use Feynman diagrams to impose conventional 
renormalization schemes such as on-shell or MS. 

It is tempting to perform a (singular) gauge transformation 



H = U^HU 



with 



U 



s e 



-mip 



-IC 



S C 



-IC 
imp 



Pl + Pr 



(11) 



such that 



ns 



(/« 



G 



-a ■ <p 
& ■ ip 



a • (p 
—a ■ p 



I G + m(f H -l) 



1 



(12) 



However, this only shifts the problem to the origin. As limp^o Ig = the - term induces a 
potential in the corresponding second order Schrodinger type scattering problem. This is outside 
the standard regime in which analyticity of the Jost function is established (25| . We comment on 
the use of H mt at the end of section V. 



1 We utilize the standard representation for the Dirac-matrices and perform a global chiral rotation U 
\ 0/ ^ >h ^ >R ^° reslore lne typical form of the Dirac equation. 
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B. Unwinding the String 

To set up the scattering problem we must find a way to treat eq. (fTUj) , A possible solution 
to that problem can be envisaged immediately. If the angle £1 were not constant, but a radial 
function going from its value at the origin (which defines the configuration of interest) to zero at 
some large distance po, the scattering problem would be treatable without altering the magnitude 
of the Higgs field. This configuration corresponds to a return string on the chiral circle that has 
less severe ultra-violet singularities, since the C3 and C4 counterterms in eq. ([5]) vanish on the chiral 
circle. For convenience we parameterize sin£i rather than £1 itself, 



1 - tanh w £ ^ 

■m = 1 + t l (m :; y - d *) = vr^w- (is) 

Of course, the chiral circle condition continues to hold. This addition does not alter the form of 
the Dirac Hamiltonian, eq. (jl(J|) . However, it must be taken into account that the isospin matrix 
Iq now has radial dependence. The combination of eqs. ([7]) and (|13p then describes the physical 
string for p <C po and a non-interacting theory for /)> p$. To accomplish the separation of scales 
indicated in eq. (P) we require po 3> wg,h- With these prerequisites, the Hamiltonian, eq. ©, has 
a well-defined scattering problem that we will now describe. 



C. Scattering off the Unwound String 

As a first step we introduce grand-spin type states that couple spin and isospin, to account for 
the angular dependence. For fixed angular momentum £ there are four of them, 



(<p;SI\£ + 



(<p;SI\1 



c 



i(£+n)<p 



(<P;SI\£ + 



(<p;SI\£ 



-1 e 



(14) 



where S and I refer to the spin and isospin subspaces, respectively. These grand-spin states serve 
to construct the Dirac spinors in coordinate space, 



(p\ + 



(p\ 



9x(p)\t-+) 



<pI + 



h{p)\£ + 

92(p)\£~ 

U(j>)\e- 

9i(p)\i + 



(15) 



where we have suppressed the angular momentum index of the radial functions because the Dirac 
equation is diagonal in this quantum number. We combine these eight radial functions into two 
vectors 



/ 



//i(p)\ 
/ 2 (P) 
MP) 

\U(P)J 



and 



fgi(p)\ 

92{p) 
93{p) 

\gi{p)J 



(16) 



to write the Dirac equation as a set of eight coupled first order linear differential equations in the 
matrix form 



(E-m)f = V uu f + (D U + V ud ) g 
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(E + m) g = (D d + V du ) f + V M g. (17) 
The derivative operators are fully contained in the diagonal matrices 

L> u = diag dp H , d p H , -d p H , -d p + - 

V p p p p 

£> d = diag -d p + , -dp + -, dp + , d p + . (18 

V P P P P / 

We will give the explicit form of the real 4x4 matrices Vi in terms of the radial functions when 
we set up the Born series for the scattering data. Here it suffices to note that these matrices are 
real and vanish at spatial infinity, so the asymptotic solutions are cylindrical Bessel and Hankel 
functions. In particular, the Hankel functions 

u u = diag^jM^^lM^i+l+^M^SlM) 

n d = diag^^^M^^^M^^UM^^lM) (19) 

that parameterize the outgoing asymptotic fields with (radial) momentum k can be used to set up 
the scattering problem via the matrix generalization 

/ — > T-U u and g — > kQ -U d . (20) 



Note that k = sgn(E') y %^ is well defined for either sign of the energy eigenvalue since \E\ > m 
for the scattering solution. Using the dispersion relation for real momenta E 2 = k 2 + m 2 , we may 
also write k = E + m = E ~^ n ■ The boundary conditions for the 4x4 complex matrices T and Q 
are simply limp-^ T = 1 and limp^oo Q = 1 so the various columns of the above products refer 
to outgoing waves in different channels. Demanding finally that the scattering wave-functions are 
regular at the origin defines the scattering matrix. It can be obtained from either the upper or 
lower components 

S = - lim %- 1 ■ J- 1 -T*-H* 

s = - lim n d x -g~ l -g* -n d . (21) 

This construction ensures that the physical scattering solution, 2 

= r ■ n* u + g ■ h u ■ s , * w = g* ■ n* u + t ■ n u ■ s . (22) 

is regular at the origin, p — > 0. An important numerical check (besides unitarity) is that both 
equations yield identical scattering matrices. We write the Dirac equation for the matrices F and 
g in a form that will later simplify the Born series, 

d p J- = [Mff + O d ] -F + F-Mf f + [M fg + kC} -g-Z d 

d p g = [M gg + o u ] ■ g + g ■ M { p g + [M gf -kc]-F-z u . (23) 



The matrices 



C = diag(-l, -1,1,1) 



Since the matrices Vi are real, T* ■ Hu and Q* ■ H* d also solve the Dirac equation. 
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O u = idiag (-{£ + n + !),-(£ + l),£ + n,> 



O d = ±diag (£ + n,£,-(£ + n + !),-(£ + 1)) , 



<l +1 (M H&ikp) H£ n (kp) Hf\kp) 



r(i) 



(l) 



M {r) 
J l gg 



diag 



< } „(M H?>(kp) H^ n+1 (k P ) H^(kp) 



(i) 



r(i) 



fcC • Z u - O u and 



(r) 
// 



-kC -Zd-Od, 



(24) 



are purely kinematic, while the background fields are contained in the real-valued 4x4 matrices 



M 



M 



<)<) 



ff 



a G R+ apR- \ 
-api?_ —a G R + J 

-a G R+ apR- 
—apR— a G R+ 



M » = L K 



M 



fg 



where 



Ra 



s c 
c —s 



R- 



a H C + a G 
a H C - a G 



-1 

1 





R + 



R+ 



-R, 


-R, 




The profile functions enter via 

olh = m (cfu — 1) , otp = msfn and a G 



ns 



fa- 



(25) 



(26) 



(27) 



Note that the inclusion of the return string makes i2+ a function of p, while R- remains constant 
in space. 

The Born series is a straightforward but tedious expansion in the matrices Alj, which are linear 
in the background fields. For example, at first order we have to solve 



d p T x = O d ■ Ji + Ji • Mf f + kC ■ Si ■ Z d + M ff + M fg ■ Z d 
d P Qi = O u -Q x +Q x - M$ -kC-F 1 -Z u + M gg + M gf ■ Z u , 



(28) 



with the asymptotic conditions lirnp_>oc F\ = and lim^oo Q\ = 0. The first order correction 
to the scattering matrix then becomes S\ = lim^o^u 1 [-^"l — J^f] H* u . The solutions F\ and Q\ 
are the sources for the second order terms Fi and Q2, yielding the second Born approximation 
52 = lim p _ s> o%M 1 [F\ (F* — F\) + Ti — F%\ H-u- This iteration in the background fields can be 
repeated to any desired order, yielding the Born series for the scattering matrix S = t+S\+S2+- ■ ■ ■ 



D. Jost Function 

Ultimately we need to compute the phase of det(<S), which according to eq. ([2Tj) is related to the 
phase of F = F^k) = det (lim^o F) , or equivalently G = Gi(k) = det(lim p _>.o Q)- Furthermore, 
we need to sum over angular momentum £ and integrate over radial momentum k after subtracting 
sufficiently many terms of the Born series. This procedure is not only numerically cumbersome 
because these are oscillating functions of k, but it may even cause erroneous results [13] since 
these sums/integrals are not absolutely convergent. These obstacles are avoided by analytically 
continuing to imaginary momenta t = ik and performing the integrals along the branch cut t > m. 
The analytic continuation for the Dirac equation is conceptually different from the well-studied 
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Schrodinger case because E = ±Vk 2 + m 2 causes the complex momentum plane to have two 
sheets. So on the real axis we have to pick one sign, continue to complex momenta and compute 
the Jost function on the imaginary axis. This procedure must then be repeated for the other sign 
and then all discontinuities must be collected at the end. In the present problem we are fortunate 
because the solutions to the Dirac equation exhibit charge conjugation symmetry along the real 
axis. Therefore det(«S) does not change under E — > —E and there is no additional discontinuity in 
the Jost function. Moreover, the Jost function is real on the imaginary axis, as in the Schrodinger 
problem. However, the way this comes about in the string problem requires us to be careful when 
constructing the Jost function for complex momenta. 

We describe the case with E = y/k 2 + m 2 . The analytic continuation concerns the Hankel 
functions, which turn into modified Bessel functions: Z u — > Y u and — > Y d , with 

_ / K e+n (tp) K e (tp) K e+n+1 (tp) K e+1 (tp) \ _ _ x 
Y U - diag y Ri+n+i{tp) , Kf+i{tp) , Ke+n(tp) , mp) J - W . { M) 

Furthermore we have the change of the kinematic coefficient k — > z K , with 

m + iy/t 2 - m 2 , nn . 
Zk = , (30) 

a pure phase. The differential equations, which we solve subject to the boundary conditions that 
T and Q approach unity at spatial infinity, read 

d p T = \M ff + O d ] -T + F-Mf f + \M fg -tC] -Q-Y A 

d P G = [M gg + O u ]-G + g- + [M gf + tC]-F-Y u . (31) 

For simplicity we have omitted the momentum arguments in the radial wave-functions J- and Q. 
It should be apparent from the context that these are distinct from the solutions to eqs. (|23p . The 
coefficient matrices are slightly modified: 



M g r) g = -tC -Y u -O u Mf f =tC-Y d -O d 



M 



of 



a H C + a G ( Q 



M 



fg 



-R^ 
R+ 



(32) 



while Mgg and Mff are as on the real axis. Note, that in contrast to the Schrodinger problem, 
the differential equations do not become real on the imaginary axis. Rather, charge conjugation 
E —> —E induces complex conjugation. It is therefore not surprising that the naive continuation to 
limpet) -T 7 is n °t real. Instead, we find numerically that F = G* . Interestingly enough, for a given 
value of t, the imaginary part is the same for all angular momenta. The origin for this imaginary 
part lies in the subtle definition of the Jost function via the Wronskian between the Jost solution, 
i.e. T or Q, and the regular solution that satisfies momentum independent boundary conditions 
at the origin. The latter are required to ensure the regular solution to be an analytic function. 
Analyticity of the Jost solution is guaranteed by the non-singular behavior of the interaction 
potentials in the presence of the return string. At the origin, the Higgs field differs from its 
vacuum expectation value (it actually vanishes). In the Dirac case this modifies the relative weight 
of the upper and lower components. On the real axis this weight is given by k in eq. (|20p . More 
precisely, at the origin the non-diagonal elements of the matrices in eq. (|32p vanish and the eight 
differential equations decouple with respect to the index on the radial functions in eq. (I16p . i.e. 
spin and weak isospin. For real momenta a typical solution in the vicinity of p = then looks 
like 



fk\ l ( V E + mc fH(0) Ji{qp) 
Q ' WE-mcfH(0)Ji+i(qp)j 



94 ' 



(33) 
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with q = \J E 2 — (mcfniO)) 2 and similar dependencies for the other six radial functions. The 
square-root coefficients cause the proper definition of the Jost function, v(t), to be 

exp [u(t)} = ( T ~7\n^ 2 lim deb (?) = ( , T + l ? f rt ) 2 lim det (^) ( 34 ) 
\T-imcf H (0)J p^o \t + imcf H {0) J p->o 

with t = Vt 2 — m 2 . The power of two occurs because we compute the determinant of a 4 x 4 
matrix. Note that this redefinition not only cancels the imaginary parts, but also modifies the real 
part. Furthermore it cancels the logarithmic singularity in In linip_j,o det ( J-) observed numerically 
at t ~ m. Since fjj is part of the interaction, this correction factor also undergoes expansion in 
the framework of the Born series. Otherwise, i.e. for det (J-"), the Born series is constructed as for 
real momenta by iterating the differential eq. (|3ip in Mi. 

The resulting Jost function is a continuous function in the upper complex momentum plane 
and the branch cuts in the Dirac equation do not carry over to v. This is a consequence of charge 
conjugation invariance in the present model. 



III. PHASE SHIFT APPROACH IN D=3+l 

In this section we collect all the pieces needed to compute the vacuum polarization energy for 
the configuration in which the physical and the return strings are combined. 



A. Interface Formalism 

The interface formalism addresses the problem of computing the vacuum polarization energy 
for configurations that are translationally invariant in a subset of coordinates. As laid out in 
ref . [3] , this approach only requires scattering data from the lower dimensional space in which the 
configuration is non-trivial. However, we still have to adapt the formalism to the present string 
problem and, in particular, motivate the above recipe for computing the logarithmic Jost function 
v(t) on the imaginary momentum axis. 

Our starting point is the interface formula for the vacuum polarization energy per unit length 
of a configuration that is translationally invariant in one coordinate, 3 



AE 



(AT) 



dk 



7T 



(k 2 + m 2 )ln 



k + m z 



I- 1 



A. 

dk 



Mk)] 



N 



(35) 



A few remarks on the notation are in order. The most important input is the phase shift 8i = 
lndet(<S) of the scattering matrix in eq. (|21|) . Using the factorization property of the determinant 
and the fact that Y v (z) 3> J u { z ) for z — > we may write 



1. 



5i(k) = ^lndet lim Tt(p, k) 1 T^(p, k) = -In det lim Ge(p, k) l Q}{p, k) . 



(36) 



The square bracket with index iV in eq. (|35p indicates that the first N terms of the Born series have 
been subtracted. This is important both for the integral to converge and also to ensure that the 



In comparison with eq.(6) of ref. [191 ] a factor two emerged from the two signs of the single particle energies. 
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semi-circle at infinity does not contribute when performing the integral by contour integration in 
the complex momentum plane. Later we will add back this contribution in terms of renormalized 
Feynman diagrams. This is indeed the heart of the phase shift approach, as it allows the imple- 
mentation of standard (perturbative) renormalization conditions. In case of the string background 
the degeneracy factor in eq. ([35]) is 



D, 



1 , £ = -n 

2 , otherwise , 



(37) 



where n = 1 is the winding number of the scalar field introduced in the string background eq. ([6]). 
Finally eji denote the bound state energies in the £ th partial wave with ejj < m. The renormal- 



ization scale p has no effect as a consequence of sum rules for scattering data [26(] . We nevertheless 
require it for dimensional reasons and for simplicity take p = m. 

To set up the contour integral we first remark that [cf. eq. ([24"]) ] Z^—k) = —Z^(k). This 
implies that 4 Fe(p,—k) = F^(p,k) since the boundary condition for F at p — > oo is independent 
from k. (Identical relations hold, of course, for Q.) We therefore write for the momentum integral 
in eq. <^ (with E = \/k 2 + 



00 dk 

7T 



2 fE 2 
E 2 ln 



dk 



Mk)] 



N 



00 dk_ 

-oo 2 ™ 

00 dk_ 

-nn 27TZ 



F 2 

f"\n ( =^ 

F 2 
m z 



— [\ndetT £ (0,k)] N 



dk 



In 



E + m 



E + mcf H (0) 



+ ]ndetJ r t (0,k) 



(38) 



N 



where the argument p = is understood as the limit p — > 0. The additional term introduced 
in the second equation adds a piece to the integrand that is odd in k and therfore does not 
change the integral. This modification has introduced the Jost function, since is precisely 
the analytic continuation of the quantity in square brackets. We have already argued that v{t) 
is an analytic function in the upper half momentum plane. Hence when closing the contour, we 
only need to consider the discontinuity of the logarithm along the cut t > m, which is 2ni. In 
addition we need to add the contributions from poles to the contour integral. As in the non- 
relativistic case, the existence of any bound state is related to a simple root of the Jost function 

at k = iKj i = i\ m 2 — (e, ^) 2 . Therefore 



tr 



F- l —F 
dk 



1 



k — iKj t £ 



(39) 



in the vicinity of k = iKj t e. The contribution of the resulting pole cancels against the explicit 
bound state contributions in eq. f)35[) . (This is a general property of the phase shift approach.) As 
a result, eq. ([35]) can be recast to 



AE 



(AT) 



— J dtt^De [u e (t)] 



N 



(40) 



4 This also determines k = 



E + m 



as the genuine starting equation for the analytic continuation in k. 



12 



Here we have integrated by parts and also used the important property that the integral along 
the imaginary axis can be interchanged with the angular momentum sum 131] . A final change of 
variables t — > r = Vt 2 — m 2 yields 



AE 



(N) 



^1 drrY.D, 



2vr 



^l(v t 2 + m 2 ) 



N 



(41) 



Equation (j41 j) is our master formula for the vacuum polarization energy of the string. However, 
we still have to make sense of the N subtractions. 



B. Feynman Diagrams in MS 



Having subtracted the first N terms of the Born series to the vacuum polarization energy, we 
must add them back in form of Feynman diagrams. We then combine these diagrams with the 
counterterms, eq. (0), to renormalize the theory. 

We generate the Feynman diagrams by expanding the effective action 



with the interaction 



and the isospin operators 



A = A - i Trln jl + (ip- m)" 1 #/} 



Hi = L^P L + h + ipj 5 



L = 0, L = 2aG<flG(f) 



h 



-awl, P=~a:plp(f) 



(42) 



(43) 



(44) 



The isospin matrices Iq and Ip are defined in eq. Q . The first two orders of the expansion contain 
quadratic and sub-dominant logarithmic ultra-violet divergences. The third and fourth order 
terms are merely logarithmically divergent. The latter divergence can be extracted by omitting 
all external momenta in comparison with the loop momentum. These are the local contributions 
to the Feynman diagrams. To completely capture the divergence structure we therefore need to 
compute the first and second order completely (to control the quadratic divergences), but we may 
restrict ourselves to the local contribution in the case of the third and fourth order terms. This 
simplifies the computation considerably. In dimensional regularization the divergent contributions 
to the action from all four orders are (where L ■ L 

d D l 



L^, etc.) 



^(div) 



Xtrji -{d a L p -i[L a ,Lp}) 2 



(h 



^2 , 2 
m) + p 



L-L 



{dhf - (dp) 2 +2{h-m)L- dp-iL- [dp,p] 



+2m 2 



(h — m) 2 + p 2 



m 



+ 



{h 



m) 2 + p 2 



m 



(45) 



The counterterms in eq. ([5]) are identical to the integrand of the spatial integral. 5 Hence these 
divergences are removed by writing the counterterm coefficient as (s = 1, . . . , 4) 

*-D f d D l 



(27T 



i-D 



(l 2 -l + ie) +c s 



(46) 



For example, tr 7 (D M $) t (L>"$) = ^ tr 7 { {dhf + (dp) 2 + [(h — m) 2 + p 2 ] L ■ L — 2 (h — m) L ■ dp + iL ■ [dp,p] }. 
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where the c s are the finite parts that are fixed by the renormalization conditions. 
The remaining ultra-violet finite contribution from the first and second order is 



i f d 4 k r 1 . r 



dx In 



(2vr) 4 

. k 2 

1 - x(l - x)-^ 



m 



6 



tr/ [h(k)h(-k) +p(k)p(-k)\ 



8?r 2 7 (2vr) 4 7 

xtr 7 |[m 2 - x(l - x)k 2 } [\L(k) ■ L{-k) - 3 (h(k)h(-k) + p(&)p(-*0)] 
+x(l - x) [k ■ L{k) k ■ L{-k) - \k 2 L{k) ■ L(-k)] 
+2m 2 p(k)p(—k) + imk ■ L(k)p(—k) > , 



(47) 



where the Fourier transforms have been indicated by the arguments of the profile functions. Observe 
that the m 3 and m 2 terms in the first two integrals combine to be proportional to an integral over 



tr 7 [$t$ 



v 



Hence they cancel exactly against the C3~type counterterm if we set its finite 



part to C3 



1 m* 



fm 2 



i -j , j 4t^' This defines the no-tadpole renormalization scheme, as the vacuum 
polarization energy is at least quadratic in the external fields. For the remaining counterterms, we 
discard all finite pieces, 



ci = c 2 = c 4 = . 



(48) 



This defines the MS renormalization scheme that we impose further on. From eq. (|47p we then 
determine the Feynman diagram contribution to the vacuum polarization energy up to quadratic 
order in the profiles to be 



AE- 



FD 



kdk k 1 



4tt 3 



{hi + p 2 n ) + Am 2 hp 2 n + 2nikha { + ) p n + k 2 I 2 



a 



(+) 



(-) 



(m 2 h + k 2 I 2 ) 



6/Iq + 6pi + ( « 



(+) 



+ a 



+ K)' 



(49) 



with the parameter integrals (77 = k/m) 

/•l 2 

h = dxki \l + x(l - xW] = - Jl + 7] 2 arsinh U) - 2, 
Jo V 

^ = f Q dx < 1 ~ x ) ln t 1 + < l ~ X W] = ^3^3 [V 2 ~ 2] arsinh (2) + JL _ A , 



(50) 



and the Fourier transforms that contain the background profiles, 



h {k) 



p dp a H (p) Jo (kp) , Pn(k) 



pdpa P (p)J n (kp) , 



«f W = / pdpa G (p)c(p) [J n -x{kp) ± J n+ x(kp)] 
Jo 



a s (k) 



pdpa G (p)s(p)Ji(kp) . 



(51) 



We still have to work out the third and fourth order pieces, for which we want to avoid computing 
the full Feynman diagrams. Essentially we only need to compensate the logarithmic divergence 

4-D r d D l . 2 . ,-2l 



A 



(div) 



3.1 



7TCF TL 



m 



(27T) 



(52) 
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where T and L are the (infinite) lengths of the time and z-axis intervals. The constant of propor- 
tionality is a simple integral over the profile functions, 

cf = J pdp | (o? H + a 2 p) (ajj + «p + 4ma#) + 4 (a 2 H + ap + 2ma.Fi) o? G 

An ~\ 

ac®. p (sap + can) j- (53) 

This result implies the limit, cf. eq. (f4T)|) . 



^D t t[u t (t)] 2 ^ -f , (54) 

I 

which is well suited to verify the accuracy of the numerical results. 

C. Fake Boson 

As already mentioned, dimensional analysis of the Dirac problem suggests that N = 4 sub- 
tractions are necessary in eqs. f)35[) and (|4ip . This is unfortunate since the Feynman diagrams for 
N = 3 and N = 4 are complicated higher dimensional integrals (in Fourier space and over Feyn- 
man parameters). Also the corresponding Born orders are numerically costly when integrating the 
Dirac equation. As discussed above, we would like to employ an alternative method to deal with 
the logarithmic ultra-violet divergences that emerge at these orders. Here we will briefly explain 
one possibility. An important prerequisite for this, however, is the above analysis on the imaginary 



axis that allows us to perform the angular momentum sum prior to the momentum integral 13]. 

Essentially we do not want to subtract the third and fourth order terms of the Born series in 
eq. (|4"T1) . but some other quantity that satisfies the following three conditions, 

1) it must exhibit the same analytic properties as the N = 3, 4 terms of the Born series, 

2) it must cancel the logarithmic divergence at large t: asymptotically it must behave like ^f-, 

3) we must be able to add it back as a sum of Feynman diagrams. 

The perfect candidate that automatically satisfies conditions 1) and 3) is the second order contri- 
bution to the vacuum polarization energy of a background potential V coupled to a fluctuating 
boson field. To ensure the second condition, we merely need to fix the strength of the background 
potential such that the divergence of the corresponding Feynman diagram matches eq. ([52)) . For 
definiteness, we take an exponential potential 

V{p) =m 2 -^e- 2p / p0 . (55) 
Po 

that also is translationally invariant along the z-axis. We choose the scale to be set by the 
return string because it determines the regions of momenta and angular momenta that dominate 
the integral and sum in eq. (|4ip . We call Vi(t) the Jost function of this boson problem on the 

(2) 

imaginary axis and v\ (t) its second Born approximation and define 



■ max 



u(t) = lim V D e [u e (t)] 2 + — lim V (2 - fc, ) vf (t) . (56) 



C-B f milx ->oo „ 
=— n £=0 
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Here 



■a-IjT^VW-?^ (57) 



is the boson analogue to eq. (j53f) . We have taken care to allow the (numerical) cut-offs for the 
angular momentum sums to be different in the fermion and boson cases. We merely have to ensure 
that, at a given t, either sum has converged to a sufficient accuracy. 
Now we can perform the momentum integral and compute 

i roo 

AE 5 = / dTTv{yJr 2 + m 2 ) (58) 

27r Jo 

as the phase shift contribution to the vacuum polarization energy. 

Finally, we find the total vacuum polarization energy of the unwound string to be the sum of 
three terms 

AE = AE S + AE FD + AE B . (59) 
The last term originates from the finite part of the second order boson contribution, 

. cf f°° kdk , T o . . 



CB Jo 16-7T 

with the Fourier transform 

r°° R — k 2 n 2 

Vb(fe) = / pdpV(p)Ukp) = m 2 pl— ^. (61) 

J ° [4 + k 2 p 2 ] * 

Two remarks on the fake boson method are in order. First, it can be used with any renormalization 
prescription because the counterterms are determined from the two leading orders as a consequence 
of gauge invariance. Second, it is important to have identical masses for fake boson loop integral 
as in the coefficient of V 2 in eq. (I59p . We note that this fake boson procedure has already been 



tested in similar models, including a consistent treatment of its mass parameter 22]. 



IV. RETURN STRING IN THE SINGULAR GAUGE 



The techniques introduced so far allow us to compute the combined vacuum polarization energy 
of the physical and return strings, but we still have to disentangle them. In general this is a tedious 
or even impossible task. However, for the case that the return string is restricted to the chiral circle, 
the gauge-transformed Hamiltonian, eq. (fl2|) . suggests a potentially successful procedure. The pure 
return string is characterized by fa = /h = 1 which apparently leads to H mt = 0. This would 
actually imply that the return string has zero vacuum polarization energy. However, this conclusion 
is premature since eq. (|12p has been derived for a constant angle £i. Repeating the same gauge 
transformation but allowing a radial dependence as in eq. (113[> yields 



H 



K a ■ p / H p\v-V / ^ 2c V a ' P ~ a ■ Pj 



Here 



/ _ ds (p) _ sin(^i) wo 
dp 1 + tanh(wo) po 



tanh 2 ( wo 



po 



(63) 
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As a consequence of chiral symmetry, the gauge-transformed Hamiltonian only has an induced 
vector field while the (pseudo) scalar fields have disappeared. 

In contrast to the gauge-transformed Hamiltonian of the physical string, eq. f)12[) . no singularity 
emerges as p — > and we have a well-defined scattering problem. It can be treated in complete 
analogy to what we described in previous sections using 



M 



M 



llll 



ff 



ay 



ay 











R- 


R- 








R- 



M 



<jf 



M 



fa 



Z K OL r 



z K a r 



R. 
R- 

R. 

R- 



(64) 



where a r = J- sets the order parameter of the corresponding Born series. Note that these are the 
coefficient matrices on for imaginary momenta, i.e. they are to be used in eq. (|3ip . 

The induced vector field only has a single isospin orientation, so the commutator term in the 
corresponding field tensor vanishes and only second order pieces contribute in eq. (|45p . As a result, 
the ultra-violet divergences from the third and fourth order Born terms vanish, which simplifies 
the renormalization procedure for the return string considerably. In particular, we do not need to 
apply the fake boson procedure. 

For the return string we do not need to take special measures to impose the no-tadpole condition 
because the C3~type counterterm vanishes on the chiral circle. Hence in the MS scheme the vacuum 
polarization energy of the return string becomes 



1 f°° 



(V 



t 2 + m 2 ) 



(65) 



where U£ is the Jost function associated with the Hamiltonian, eq. (I62j) . and 



AE- 



(r.s.) 
FD 



kdk r 
IvTl 



m 1 1 



a 



+ 2k 2 I 2 («(")) 2 } 



(66) 



is the finite piece of the second order Feynman diagram. It includes the Fourier transforms 



a^(k) = / pdpa r (p) [J n+1 (kp) ± J n -x(kp)} . 
Jo 



(67) 



We are now in a position to present numerical results. 



V. NUMERICAL ANALYSIS 

In this section we present results for the vacuum polarization energy 

E ps = AE — A£ (r - S -) (68) 

of the non-abelian string generated by the fluctuating fermions. The fermion mass m = vf sets 
the scale for all data presented in this section. The length parameters wg,h and po are measured 
in units of ^ , and since E ps is the vacuum polarization energy per unit length of the string, it is 
measured in m 2 . Note that the inverse width parameter (wo) associated with the return string is 
dimensionless. As mentioned earlier, all results presented here are for unit winding of the string, 
n = 1. 

The main purpose of the present study is to demonstrate, as a proof-of-principle, that the 
leading (fermion) contribution to the vacuum polarization energy of a string in a non-abelian 
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Po 


6 8 10 12 14 


£0 


0.0065 0.0024 0.0014 0.0010 0.0008 



TABLE I: Lowest bound state energy for £i = wh = w G = 2 and w = 8 as a function of the position of 
the return string, p . 

gauge theory can indeed be unambiguously computed. This has made necessary the introduction 
of a return string. We parameterized its shape by two variables, the position po an d the (inverse) 
width wo. In this study we concentrate on the variation of E ps on these two parameters, to establish 
a well-defined limit. The underlying feature of our calculations is the separation of the physical 
and return strings' vacuum energies. Before presenting the corresponding results, we like to discuss 
a necessary condition for this separation to be observed. The physical string generates a zero mode 
for £i = ^ in the i = — n channel 0], regardless of the specific shape that is parameterized by 
wh and wg- This zero mode must also be produced when the physical string is augmented by 
the return string if the separation hypothesis holds, because the bound state wave function does 
not interfere with the return string. In table U we present the lowest single particle energy 6 for 

= | as a function of the separation between the physical and the return string. As expected, 
this eigenvalue approaches zero as the separation increases. These results provide confidence in the 
separation hypothesis already at moderate values of po- 

The numerical calculation of the vacuum polarization energy is quite involved. The main effort 
concerns the phase shift part, AE$ in eq. ([59]). In addition to the criteria mentioned after eqs. (|21|) 
and (|54p we have furthermore tested these calculations of scattering data with respect to 

• charge conjugation symmetry: E -4 —E (for real momenta) 

• reflection around the winding of the string: I — >■ —i — n 

• sum rules for distinct angular momenta, i.e. agreement of integrals involving the Jost- 
function over real and imaginary momenta. 

We integrate the differential eqs. ()23p . (|3ip . their Born expansions as well as the fake boson and 
return string analogues from some large radius p max ~ 4po to p m ; n ~ with the boundary condition 
J(/) m ax, k) = 1, and identify lim p _>.o F(p, k) = T(p m ^,k). Alternatively, this identification can 
equally well be obtained from the derivative of the wave-function. Furthermore, a differential 
equation is formulated for lndetJ-"(p, k)~ 1 J r (p, k)* to avoid 2it ambiguities in the computation of 
the phase shift, 5e(k), cf. eq (|36|) . The computations for real momenta have been performed mainly 
for use in the consistency tests mentioned above. Channels that include (modified) Hankel functions 
with zero index (I = —2, —1, 0) are particularly cumbersome because in separating the regular and 
irregular solutions, we need to distinguish ln(p) from a constant at very small distances p m i n , cf. 
eq. (i22|) . That is, /9 m i n must be taken tiny to obtain the correct scattering matrix in eq. (f21~|) . On 
the real axis, the result can be checked against extracting the 5-matrix from the derivative of the 
scattering wave-function. For calculations on the imaginary axis we assume p m - m ~ 10~ 60 and 
successively carry out an extrapolation 

KPmin) = + ; ^ + T-^-T — V • • • ' ( 69 ) 
ln(p m in) ln^(p min ) 



6 We compute this energy by discretizing the spectrum of the free Dirac Hamiltonian in a large cylindrical cavity, 
computing the interaction matrix elements from those eigenstates, and diagonalizing the resulting matrix. 
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TABLE II: Numerical results for the vacuum polarization energy of the return string in the MS scheme. The 
entry A£ , ^ r s ' ) refers to the t integral in eq. (|65p. It is the analogue of eq. (|58|) for the sole return string. Scales 
are set by the fermion mass m. 

for the Jost function in these channels. We test the final result, i.e. uq, for stability against further 
changes of p m i n and also check the condition lm(fo) = 0. In other channels p m i n ~ 10~ 12 suffices 
to represent the origin. 

For the sum over angular momentum channels, we must go to very large channel numbers 
(typically several hundred). The required value ^ max to obtain convergence increases with po. 
While this is expected, since po determines the impact parameter, it limits our ability to send ^ max 
to infinity. We therefore perform an extrapolation to £ max —> oo. Using various analytical forms 
for this extrapolation we estimate an absolute error of 0.005 for AEg. Fortunately, the right hand 
side of eq. (j56[) approaches zero at moderate momenta, so that we can approximate the upper limit 
of the integral in eq. ()58[) by 10 or even less. 

The individual sums in eq. (|56p have typical maximal values of the order of tens, while they 
combine to a result at the order of a tenth or less. This behavior is typical for these type of com- 
putations because the physics information is hidden beneath the ultra-violet divergences. This 
disproportion increases with the distance of the return string. For practical purposes, this unfortu- 
nately implies a severe loss of precision of three orders of magnitude or more. However, demanding 
too high numerical precision makes the integration of the Dirac equation (and its Born expansion) 
ineffective. Fortunately, doing the momentum integral on the imaginary axis regains quite some 
efficency because we are not plagued by oscillating integrands. 

In total, it takes several CPU-days on a modern PC to generate a number for AE for a single 
set of parameters (wg,h, Ci, wo, po). The fake boson piece YLl^li m comparison, is cheap since it 
takes less than a CPU-hour even though l max ~ 3^ max for sufficient convergence. 

We first consider the return string contribution, which only depends on the parameters Wo and 
po- The results are shown in table [III We observe that AE^ T,S ^ is small in magnitude and quickly 
decreases as the position of the return string is sent to infinity. Dimensional analysis indicates that 
except for the Higgs kinetic term (the C2~type in eq. ([5])) all contributions from the return string 
vanish as po - > oo. This is the local part of the second order Feynman diagram, which is taken 
out in the MS scheme. Hence both the smallness and the large po behavior of AE( T - B -1 are well 
understood. Since for our particular parameterization, the relevant parameter is wo/po, the above 
argument also explains the decrease with wo that can be observed from table [Hi In essence, the 
contribution of the sole return string is well under control. 

In table IIIII we show the contributions to the vacuum polarization energy of the configuration 
consisting of both the physical and the return string, for various parameters that describe the 
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TABLE III: Numerical results for the vacuum polarization energy of the combination of the physical and the 
return strings in the MS scheme with wh = »g = 2 and £i = 0.47T. 
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TABLE IV: Vacuum polarization energy of the physical string for wh = 2, wq = 2 and £i = 0.47T. 

return string part of the configuration. We observe a significant cancellation between the integral 
over the subtracted Jost-function and the fermion Feynman diagram contribution. This behavior, 
of course, can be different in renormalization schemes other than MS. The fake boson Feynman 
diagram contribution is negligible in the MS scheme. 

In table |IV] we combine the results from the previous two tables to find the vacuum polarization 
energy of the physical string. Certainly the data indicate a saturation for large po, as demanded 
by the assumption that the scales for the physical and return strings separate This separation 
assumption is further supported by the observation that for increasing distance of the return 
string, its width (^) no longer affects the energy of the physical string. In this limit, the vacuum 
polarization energy of the string turns out to be less than a tenth of the fermion mass. We have 
repeated this analysis for a thin string and present the results in table |V] We find that its vacuum 
polarization energy is even smaller. 

We recall that for typical model parameters, i.e. gauge and Yukawa couplings 7 of the order 
unity or less, the classical energy of the physical string configuration, which we did not discuss here, 
is of the order 1 or 10. Also typical single particle energies of the fermions bound in the string 
background are around 0.5. Hence the vacuum polarization energy in the MS renormalization 



7 The relevant scale is set by the heaviest fermion coupling to the string, i.e. the top quark, where / ~ 175/190 ~ 1 
as the Weinberg angle is set to zero. 
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Po 


AE 


-A£( r - S -) 


AE ps 


6 


0.003 


0.107 


0.110 


8 


0.022 


0.052 


0.074 


10 


0.029 


0.030 


0.059 


12 


0.027 


0.020 


0.047 


14 


0.025 


0.013 


0.038 



TABLE V: Same as table[IV]for wh = 0.5, wq = 1 and £i = 0.47T. Here we used wq = 8. 



scale is about an order of magnitude smaller than other relevant energy scales in the model. This 
smallness of string (vortex) vacuum polarization energies has also been observed for string type 
configurations in simpler models [8] and hence is not fully unexpected. 

We have also performed numerical calculations with the singular Hamiltonian of eq. (|12p by 
regularizing the singularity of the gauge coupling: - — >■ In the second order scattering 

problem the unregularized form induces a \ divergence at small p, which eventually spoils the 
analytic properties of the scattering data. However, we did not obtain a smooth limit as we removed 
the regulator. This is not a contradiction to the existence of a well defined vacuum polarization 
energy. It just implies that two operations of (i) integrating to p — > and (ii) removing the regulator 
may not be exchanged. This can already be observed from computing the classical energy with 
such a regulator. The spatial integral, which involves the field strength tensor squared, has no 
well-defined limit as the regulator w reg is sent to zero. 

Finally, we would like to mention that we have also considered the numerical computation 
on the real axis, e.g. using eq. (I35p . to compute the vacuum polarization energy for this Dirac 
problem. As expected from simpler models, we find that a twofold subtraction suffices to render a 
convergent integral when the angular momentum sum is carried out first. This result contradicts 
the ultra-violet divergence structure of the quantum field theory. It is a mathematical artifact 
stemming from the improper exchange of sum and integrals 13j. This has been the major reason 
to analytically continue to imaginary momenta, since it allows this exchange. We remark that a 
rigorous Born subtraction, i.e. up to at least fourth order in Mi in eq. (|23p . can indeed be carried 
out on the real momentum axis. As argued before this seems technically problematic. It is the 
exchange of sums and integrals that requires absolute convergence, which is only guaranteed for 
imaginary momenta. 



VI. CONCLUSIONS 



We have presented a feasible method for the computation of the fermion contribution to the 
vacuum polarization energy of a string type configuration in a non-abelian gauge theory. If there 
are many internal degrees of freedom for the fermions, for example color, this contribution is the 
dominant quantum correction to the energy of the string. The key tool that we employed was the 
introduction of a return string configuration on the chiral circle. It was, in particular, indispensable 
to obtain a well-defined Feynman series, which is required to impose conventional renormalization 
conditions. For the present study we used the MS scheme, augmented by the no-tadpole condition. 
The generalization to other schemes, such as on-shell, does not cause further (technical) subtleties. 
The advantage of having the return string on the chiral circle is that its vacuum polarization energy 
can be reliably computed in a specific gauge. Then it suffices to compute the energy of any physical 
string configuration with respect to this return string. Furthermore, technical innovations such as 
the computation of momentum integrals via analytic continuation and the introduction of a fake 
boson field have been necessary to reduce the computational workload, making the investigation 
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feasible with appropriate expenditure of computational resources. 

Our results suggest that for future investigations of the total energy of string configurations 
it is a good approximation to omit the MS contribution to the vacuum polarization energy. This 
contribution is about an order of magnitude smaller than other typical energies in the model. This 
is in particular the case for a thin string. Such an approximation simplifies the computation of the 
energy of the string tremendously: It merely remains to find the (finite pieces of) the countert- 
erm coefficients, since the spatial integral of these terms is, by the definition of renormalizability, 
the same as for the classical energy. The remaining spatial integrals only involve gauge invariant 
combinations and thus can immediately be computed without the return string. This result will 
make the search for configurations feasible that are stabilized on purely dynamical grounds. For 
instance, such a stabilization can occur by trapping fermions along the string. Once such a config- 
uration is discovered, it remains to verify that its vacuum polarization energy in the MS scheme 
indeed is tiny, which can be achieved with the techniques presented here. We leave this study to 
a forth-coming paper. 
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